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Abstract. We present the results of nucleosynthesis calculations based on multi-dimensional (2D and 3D) hydrodynamical 

simulations of the thermonuclear burning phase in type la supemovae (hereafter SNIa). The detailed nucleosynthetic yields 

OO . of our explosion models are calculated by post-processing the ejecta, using passively advected tracer particles. The nuclear 

Cn ' reaction network employed in computing the explosive nucleosynthesis contains 383 nuclear species, ranging from neutrons, 

^•-^ protons, and a-particles to ''^Mo. Our models follow the common assumption that SNIa are the explosions of white dwarfs that 

have approached the Chandrasekhar mass (M,.;, ~ 1.39), and are disrupted by thermonuclear fusion of carbon and oxygen. But 

in contrast to ID models which adjust the burning speed to reproduce lightcurves and spectra, the thermonuclear burning model 

applied in this paper does not contain adjustable parameters. Therefore variations of the explosion energies and nucleosynthesis 

1-^ ' yields are dependent on changes of the initial conditions only. Here we discuss the nucleosynthetic yields obtained in 2D and 

t^l 3D models with two different choices of ignition conditions (centrally ignited, in which the spherical initial flame geometry 

Q . is perturbated with toroidal rings, and bubbles, in which multi-point ignition conditions are simulated), but keeping the initial 

Vh ' composition of the white dwarf unchanged. Constraints imposed on the hydrodynamical models from nucleosynthesis as well 

jy^ , as from the radial velocity distribution of the elements are discussed in detail. We show that in our simulations unburned C and 

C^ ■ O varies typically from ~40% to ~50% of the total ejected material. Some of the unburned material remains between the flame 

plumes and is concentrated in low velocity regions at the end of the simulations. This effect is more pronounced in 2D than in 

3D and in models with a small number of (large) ignition spots. The main differences between all our models and standard ID 

computations are, besides the higher mass fraction of unburned C and O, the C/0 ratio (in our case is typically a factor of 2.5 

?— j I higher than in ID computations), and somewhat lower abundances of certain intermediate mass nuclei such as S, CI, Ar, K, and 

Ca, and of '*Ni. We also demonstrate that the amount of '*Ni produced in the explosion is a very sensitive function of density 

and temperature. Because explosive C and O burning may produce the iron-group elements and their isotopes in rather different 

proportions one can get different '*Ni-fractions (and thus supernova luminosities) without changing the kinetic energy of the 

explosion. Finally, we show that we need the high resolution multi-point ignition (bubbles) model to burn most of the material 

in the center (demonstrating that high resolution coupled with a large number of ignition spots is crucial to get rid of unburned 

material in a pure deflagration SNIa model). 

Key words, hydrodynamics - nucleosynthesis, nuclear reactions ~ supemovae: general 

1 . Introduction II supemovae (SNII), SNIa progenitors are thought to be white 

dwarfs (WDs) in binary systems (see Whelan & Iben 1973, 

Type la supemovae (SNIa) are known to be stellar explosions ^nd HiUebrandt & Niemeyer 2000 for a more recent review), 

with no signs of hydrogen and helium in their spectra, but inter- i^ the canonical model the WD, expected to consist mainly of 

mediate mass elements such as Si, S, Ca and Mg near the max- carbon and oxygen, approaches the Chandrasekhar mass (Meh) 

imum of their light curves, and many Fe lines at later times, through a not yet known mechanism, presumably accretion 

In contrast to massive stars which are the progenitors of Type f^^^ ^ companion star, and is then dismpted by a thermonu- 

clear explosion. The declining light is powered by the radioac- 

Send offprint requests to: C. Travaglio live decay of ^''Ni. A strong argument in favor of this scenario 
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is given by the fact that these explosion models fit quite well 
the observed light curves and spectra (Leibundgut 2001). 

Despite the consistency of this general framework with 
observations the detailed theory of how SNIa evolve and ex- 
plode is still subject of considerable efforts. Over the last three 
decades, one-dimensional spherically symmetric models have 
been used to study the various channels that may give rise 
to a successful SN la in terms of the predicted spectra, light 
curves, and nucleosynthesis. Much of this work was centered 
on the Mch scenario wherein a C+O white dwarf accretes 
H or He from a binary companion (Nomoto, Thielemann, & 
Yokoi 1984) and ignites explosive carbon burning just before 
it reaches a critical mass of Mch ~ 1-39 M©. The subsequent 
explosion produces enough ^^Ni (~ 0.6 Mq) and intermediate 
mass elements to reproduce "normal" SN la lightcurves and 
spectra, provided that the amount of Ch-O burned at any given 
density is suitably chosen. In ID models this can be achieved 
by parameterizing the thermonuclear flame speed and, if de- 
sired, the density at which a transition to supersonic burning 
(detonation) occurs (Khokhlov et al. 1999; Niemeyer 1999). 
Moreover, some mixing of processed matter had to be assumed 
in order to fit the observed spectra. Alternative scenarios, in- 
cluding sub-Mch explosions and merging white dwarfs (double 
degenerates), have met with mixed success (see e.g. Arnett & 
Livne 1994). 

More recently it has become possible to perform multi- 
dimensional 2D (Livne 1993; Reinecke, Hillebrandt, & 
Niemeyer 1999; Lisewski et al. 2000) and 3D (Reinecke, 
Hillebrandt, & Niemeyer 2002a,b; Gamezo et al. 2003) sim- 
ulations of an exploding Mch-white dwarf. The principal diffi- 
culty in these models is the fact that the hydrodynamically un- 
stable and turbulent nuclear flame front develops structures on 
much smaller length scales than can numerically be resolved. 
However, this problem can be overcome by "large eddy sim- 
ulations", i.e. by employing subgrid-scale models for the un- 
resolved scales that provide a guess of the effective turbulent 
flame speed on the scale of the computational grid (Niemeyer 
& Hillebrandt 1995a; Ropke, Niemeyer, & Hillebrandt 2003). 
In this flame model, which is well justified in the thin-flame 
regime and is tested in experiments with premixed turbu- 
lent chemical flames, we do not need a detaUed prescrip- 
tion of the nuclear reactions. Instead, the fuel consumption 
rate is roughly propartional to the surface area of the flame 
front and its normal (turbulent) velocity. 

Despite of the need of more detailed studies of such 
subgrid-scale models, it is important to stress that the multi- 
dimensional simulations reach a qualitatively different level of 
predictive power than ID models. In particular, the amount 
of material burned at a given density can not longer be fine- 
tuned but is determined by the fluid motions on the resolved 
scales and a particular choice of the subgrid model (Reinecke et 
al. 2002a). Therefore, once the flame model has been fixed nu- 
merical simulations of the thermonuclear explosion of a given 
white dwarf can be done by just choosing the ignition condi- 
tions, including the chemical composition of the WD, the only 
remaining (physical) parameter. 

The undeniable influence of SNIa explosions on, among 
others, the chemical evolution of galaxies makes the quest 



for solid theoretical models and nucleosynthetic yields an ur- 
gent task. Guided by decades of modeling and nucleosyn- 
thesis calculations in spherically symmetric models (the pro- 
toype being the W7 model of Nomoto et at. 1984, Iwamoto 
et al. 1999, Brachwitz et al. 2000), we have begun analyzing 
the detailed nucleosynthetic yields of our explosion models by 
post-processing the ejecta. This has been performed adding a 
"lagrangian component" to our Eulerian scheme in the form of 
tracer particles passively advected with the flow in the course 
of the Eulerian calculation. Therefore we record their T and 
p history by interpolating the corresponding quantities from 
the underlying Eulerian grid. A similar method of tracer par- 
ticles in an Eulerian code to calculate the nucleosynthesis has 
been adopted in a previous study of multi-dimensional nucle- 
osynthesis in core collapse SNe by Nagataki et al. (1997), and 
more recently in calculations for very massive stars (Maeda et 
al. 2002), for core collapse SNe (Travaglio et al. 2004), and for 
Type la SNe first preliminary results have been discussed by 
Niemeyer et al. (2003). 

In this paper we present the nuclear yields resulting from 
several of our multi-dimensional supernova simulations, and 
we compare them to the standard W7 (Iwamoto et al. 1999, 
Brachwitz et al. 2000, Thielemann et al. 2003) results. In 
Section 2 we summarize the 2D and 3D SNIa calculations, dis- 
cussing different mode of ignition as well as grid resolution of 
the hydrodynamic code. In Section 3 we describe our method 
to perform nucleosynthesis calculations and the nucleosynthe- 
sis network adopted for this work. Finally, in Section 4 we 
present and discuss our nucleosynthesis results. In a first step 
we have performed resolution studies in 2D consisting of differ- 
ent methods how to distribute the tracer particles, the number of 
particles used, and grid resolution of the hydrodynamic code. 
Although 2D simulations cannot be considered to be realistic, 
as was discussed by Reinecke et al. (2002a), they can serve 
to guide the more elaborate 3D models. We then discuss the 
nucleosynthesis resulting from three 3D models, a centrally ig- 
nited model and two models with a few and many off-center ig- 
nition spots, respectively. It will be shown that the more realis- 
tic ignition conditions (central ignition or many ignition spots) 
also predict nucleosynthesis yields closer to the ones observed 
in typical SN la's. 

Concerning nucleosynthesis we will in particular analyze 
the range and distribution of ^*'Ni masses we are able to produce 
with our present models, and the sensitivity of the amount and 
velocity distribution of unburned material ('^C, '^O, ^^Ne) to 
the ignition conditions of the explosions which are still a major 
uncertainty of SN la models. 



2. Recent multi-dimensional SNIa calculations 

We have carried out numerical simulations in 2D and 3D, for 
several diff'erent ignition conditions, and for diff'erent numeri- 
cal resolution. Details of these models are given in a series of 
papers (Reinecke et al. 1999, 2002a, 2002b). A detailed dis- 
cussion of these models is not the aim of this paper, therefore 
only a summary of the results essential for our nucleosynthesis 
calculations will be repeated in this section. 
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As long as the evolution of the white dwarf before the 
thermonuclear runaway remains largely unexplored (see recent 
work by Woosley, Wunsch, & Kuhlen 2004 and references 
therein), only very crude constraints can be put on the flame 
geometry at the onset of the burning. It appears likely that the 
deflagration sets in at the surface of quietly burning "hot bub- 
bles". Nevertheless very little is known about the number, size 
and radial distribution of these hot spots. This is a consequence 
of the complicated physical processes taking place in the white 
dwarf's core during the convective smoldering phase prior to 
ignition lasting for ~1000 years. The long time scales com- 
bined with the relatively slow convective motions make nu- 
merical simulations of this phase a daunting task which has not 
been undertaken in its full complexity so far. Theoretical con- 
siderations and simplified simulations carried out by Garcia- 
Senz & Woosley (1995) suggest that fast burning starts on the 
surface of many small bubbles (r < 5 km), within 100 km 
of the star's center. Central ignition is another possible sce- 
nario that has been investigated during the last years using mul- 
tidimensional calculations (Niemeyer & Hillebrandt 1995b; 
Khokhlov 1995). 

In this work we present four models: c3J2d-512, a 2D 
model with central ignition, and grid size of 512^; c3^d-256, 
a 3D model with central ignition, and grid size of 256^; 
b5_3d^56, a 3D model with ignition in 5 bubbles, and grid 
size of 256^; finally, b303dJ768, a 3D model with ignition in 
30 bubbles, and grid size of 768^. This last one is the model 
with the highest resolution possible to evolve with the com- 
puter resources available to us, therefore we will consider it as 
the 'standard' model for this paper It achieves a central reso- 
lution of 3.33 km, using a grid consisting of 768^ zones. In the 
simulated octant of our model b30_3d_768, 30 bubbles with a 
radius of 10 km were distributed randomly. The bubble loca- 
tions were drawn from a Gaussian probability distribution with 
a dispersion of cr = 75 km. Bubbles located more than 2.5 cr 
away from the center were rejected. In all models we started the 
simulations with a central ignition density of 2.9x10^ gr/cm^. 
The simulations have been followed up to 1.5 sec. for all the 
models, except for the b30_3dJ768. Due to a very high con- 
sumption in computer time the model b30^dJ768 was stopped 
when no further energy was released. 

Fig. 1 shows the energy release for the four models men- 
tioned above. The curves are nearly identical during the first 
~0.5 sec of the simulation. Owing to the small volume of the 
bubbles, the initial hydrostatic equilibrium is only slightly dis- 
turbed. During the first stages the energy release is therefore 
lower than in previous simulations. Only after the total flame 
surface has grown considerably (mostly by deformation of the 
bubbles), vigorous burning sets in. In the late explosion phase 
(after about 0.5 sec) the total energy differs for the four sim- 
ulations, and increases moderately with increasing resolution 
(see Reinecke et al. 2002a,b for more detailed discussion). We 
also note that the 3D centrally ignited and the five-bubble mod- 
els are remarkably similar, even if the centrally ignited has a 
relatively faster burning between ~0.5 and ~1 sec. 

It must be noted that the five-bubble model is not identical 
to the model b5JdJ256 presented by Reinecke et al. (2002b): 
due to an oversight during the simulation setup the initial po- 
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Fig. 1. Total energy evolution for the two-dimensional cen- 
trally ignited explosion model {solid line), for the three- 
dimensional low-resolution centrally ignited {dashed line) ex- 
plosion model, for the three-dimensional low-resolution 5 bub- 
bles {dotted-dotted-dashed line), and for the three dimen- 
sional high-resolution 30 bubbles {dotted-dashed line) explo- 
sion model. 

sitions of the burning bubbles are not the same. As a conse- 
quence, the total energy releases of these two simulations are 
slightly different. 

The initial configuration of the front, as well as snapshots 
of the front evolution at later times are shown in Fig. 2 and 
Fig. 3 for the models c3^d-256 and b303dJ68, respectively. 

In addition to the total energy release, the mass fraction of 
unburned material in the central region of the remnant appears 
to be a good criterion for judging the validity of our simula- 
tions, because a high amount of C and O in this region would 
most likely produce a characteristic signature in the late-time 
spectra which has not yet been observed. In this respect the 
results of our earlier calculations were not very encouraging 
since the ashes rose towards the surface in large structures and 
left nearly pure fuel in the center. Using many initial bubbles, 
however, seems to alleviate this problem insofar as the statisti- 
cal isotropy of the initial flame at least delays the development 
of large-scale turnover motions. As a consequence, C and O 
is lower than 20% in the central 0.2 Mq after 0.9 s for model 
b30^dJ768. For this last model we get ~40% of the total mass 
stays unburned (we define unburned the material with T < 
1.5x10^ K). We will discuss this point in more detail when we 
will present our nucleosynthesis results. 

3. Nucleosynthesis in multi-dimensional SNIa 

The multidimensional SNIa simulations described in Section 2 
employed a minimal nuclear reaction network, directly in- 
cluded in the hydrodynamic code, sufficient for a good ap- 
proximation of the thermonuclear energy release and the pre- 
dicted gross chemical composition agrees well with the expec- 
tations (Reinecke et al. 2002b). It consists of five species {a- 
particles, '^C, '^O, ^'^Mg, and ^''Ni) and is intended to model 
the energy release of the thermonuclear reactions only. No re- 
action rates are employed: all material behind the flame front 
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Fig. 2. Snapshots of the front evolution for the centrally ignited model c3_3dJ256 at s, 0.2 s, 0.4 s, and 0.6 s. 



is instantaneously converted to a NSE of ^^Ni and ff-particles 
at high densities and to ^'*Mg at intermediate densities. Below 
lO^g/cm"', no burning takes place. 

Therefore we follow a minimal reaction network di- 
rectly in the hydrodynamic simulations, and a much more 
extended network in a post-process step. Here we present the 
results of the more detailed study of the nuclear abundances in 
the ejecta obtained by post-processing the output of the four 
hydrodynamic models discussed above. 

Since the multidimensional hydrodynamics scheme used 
in modeling the explosions is an Eulerian one (i.e. the grid 
does not move with the fluid), in order to record temperature 
and density evolution as a function of time (the necessary in- 
put for nucleosynthesis calculations) we homogeneously dis- 
tributed ~ 10000 marker particles (in 2D models) and ~20000 
marker particles (in 3D models) and followed their T and p evo- 
lution. The number of particles in the simulation was chosen in 
order to reproduce in the best way the resolution of the grid 
(see discussion below). We then calculated the nucleosynthesis 
experienced by each marker and computed the total yield as a 
sum over all the markers, after the decay of unstable isotopes. 

3.1. Tracer particles method 

In one spatial dimension it is nowadays possible to solve re- 
action networks consisting of hundreds of species online with 



the hydrodynamics (see e.g. Rauscher et al. 2002 for explosive 
nucleosynthesis calculations in core collapse SNe). However, 
it is more common to use reduced networks in order to ob- 
tain the (approximate) energy generation rate for the hydrody- 
namics and to calculate the detailed chemical composition only 
afterwards in a post-processing step. This is facilitated by the 
lagrangian nature of nearly all ID codes employed for explo- 
sive nucleosynthesis calculations. In lagrangian schemes, the 
grid moves with the fluid and therefore it is possible to record 
the evolution of the temperature and density for different fluid 
elements (i.e. lagrangian zones), which is required for the post- 
processing. In contrast, most grid-based multidimensional hy- 
dro schemes are of Eulerian type (i.e. the grid is fixed in space). 
To obtain the necessary data for the post-processing we added 
a "lagrangian component" to our Eulerian scheme in the form 
of marker particles that we passively advect with the flow in 
the course of the Eulerian calculation, recording their T and 
p history by interpolating the corresponding quantities on the 
underlying Eulerian grid. 

In the 3D simulations the star is subdivided into 27^ grid 
cells equidistant in the integrated mass M(r), azimuthal angle 
ifi and cos 0, so that each grid cell contains the same mass. 
A tracer particle was placed randomly in each of those grid 
cells, therefore the total number of tracer particles we used is 
27^ = 19683. After numerical inversion of the function M(r), 
the (M,ip,cos 9) coordinates are mapped onto the Cartesian grid. 
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Fig. 3. Snapshots of the front evolution for the floating-bubble model b30_3d_768 at s, 0. 1 s, 0. 14 s, and 0.2 s. 



For the 2D simulation, 100^ particles are distributed in r and 
cos 6 directions, using the same procedure as above. In all cases 
the simulation covers one octant, therefore to get the total white 
dwarf mass we mutiply the mass of each marker by 8 and we 
sum over all the markers. They are distributed in the way to 
have identical mass, that is therefore calculated as the ratio be- 
tween Mch and the total number of tracer particles. The mass 
of each marker is therefore between ~1Q'^Mq and ~10^'^Mq. 
The initial distribution of the markers for the b30^dJ768 case 
is shown in Fig. 4. 

Finally, we compare the r(t) values for each marker ex- 
tracted from the hydrodynamical model with the T(t) derived 
from using the internal energy from the hydro-code (which 
has contributions from the Boltzmann gas of ionized nuclei, 
the Planck-spectrum of photons, and relativistic degenerate or 
non-degenerate electrons positrons). This is done by calculat- 
ing separately the equation of state for a given maker particle's 
internal energy, density and composition, and deriving from 
that the temperature T. As chemical composition for this cal- 
culation we use a 16 isotopes network, composed by n, p, ^He, 
i^C, i^o, 20Ne, 24Mg, 28Si, 32s, ^^Ar, ^OCa, 44Ti, ^^Cr, '^Fe, 
^''Ni, ''°Zn. We found that the r(t) obtained directly from the 
hydrodynamic model is in average lower by 10% up to 20% 
for markers with high temperatures (T > 6x10' K) as com- 
pared with the temperature derived from the energy. This can be 
understood from the fact that in the hydro-code the internal en- 




Fig. 4. Radial distribution of the tracer particles in the 3D 
model at the beginning of the simulation. 



ergy density is the directly computed variable and, therefore, is 
more accuratly determined than the temperature. Consequently 
the more precise r(t) distribution is derived from the internal 
energy and the actual composition. This is what we used for 
our nucleosynthesis calculations. 
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3.2. Nucleosynthesis network 

The nuclear reaction network employed in computing our post- 
process explosive nucleosynthesis calculations contains 383 
nuclear species ranging from neutrons, protons, and a-particles 
to molibdenum. A detailed description of the code we used 
to solve the nuclear network and the reaction rate library uti- 
lized is given by Thielemann et al. (1996) and Iwamoto et 
al. (1999). Weak interaction rates applied in those calculations 
were taken from Fuller, Fowler, & Newman (1985). More re- 
cently full large-scale shell model calculations for electron cap- 
ture andjS-decays became available also forpf-shell nuclei, i.e. 
the Fe-group (Langanke & Martinez-Pinedo 2000, Martinez- 
Pinedo et al. 2000). They have already been included in pre- 
liminar calculations by Brachwitz et al. (2000) and Thielemann 
et al. (2003). We also included the new rates in the calcu- 
lations presented in this paper. As discussed below more in 
details, the nuclear reaction rates entering the thermonuclear 
modeling can play an important role. While large portions of 
the ejecta which experience maximum temperatures in excess 
of 6x10^ K follow nuclear statistical equilibrium (a chemical 
equilibrium of all strong and electromagnetic reactions), weak 
interactions occur on a longer timescale and different choice 
of Fueller et al. (1985) (as used by Iwamoto et al. 1999) or 
Langanke & Martinez-Pinedo (2000) (as used for this work, 
and by Brachwitz et al. 2000, Thielemann et al. 2003), can 
strongly affect the results. 

The initial WD composition we used consists of (mass frac- 
tion) 0.475 Mq of i^c, 0.5 Mo of ^^O, and 0.025 Mq of ^^Ne 
(in agreement with the W7 initial composition, Iwamoto et 
al. 1999). With this initial composition we typically simulate 
a solar metallicity SNIa. 

When the flame passes through the fuel, ^^C, '*'0 and ^^Ne 
are converted to ashes with different compositions depending 
on the intial T and p. We stop our nucleosynthesis calculations 
after ~1.5 sec. When the temperature in the markers dropped 
at ~1.5 10'' K the explosive nucleosynthesis is almost frozen. 
The distribution of T and p as a function of time is shown in 
Fig. 5 for the model b30^dJ768. The thick lines represent the 
upper envelope for T and p, and with the dashed lines we plot 
T and p histories in some markers taken as examples. As one 
can notice from Fig. 5, the T distribution is not yet below ~1.5 
10^ K for an the tracer particles. At T ~ 3 - 4 10^ K (i.e. the 
upper values shown in Fig. 5) we still expect some explosive 
C-burning products. Therefore we linearly extrapolate T andp 
until all the particles have T < 1.5 10^ K (that corresponds to ~ 
1 sec). 

The combination of T and p vs. time in each marker is very 
important for the nucleosynthesis results (as discussed below). 
Comparing the model shown in Fig. 5 with the T and p dis- 
tribution shown by Iwamoto et al. (1999) for their W7 model, 
we note differences that can be interesting for the nucleosyn- 
thesis calculations. First, the timescale in our models are much 
smaller (~1.5 sec) with respect to the W7 timescale (~6 sec). 
Then different combinations of T and p, in our case rather low 
T at still high p, also give us interesting differences in the nu- 
cleosynthesis calculations. 
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Fig. 5. Temperature (upper panel) and density (lower panel) 
history in the tracer particles for the b30JdJ68 model. The 
thick lines represent the upper envelope of the distribution and 
the thin dashed lines represent some of the markers randomly 
taken as examples. 



4. Discussion and results 

In this Section we present the results for nucleosynthesis cal- 
culations in the model c3J2d-512, c3_3d-256, b5^d-256, and 
our 'standard' b303d-768. We also compare them with the 
W7 calculations by Brachwitz et al. (2000) and Thielemann 
et al. (2003) (note that in Figures and Tables we better com- 
pare with Brachwitz et al. calculations instead of Iwamoto et 
al. 1999, in order to be consitent with the use of elecetron cap- 
ture rates). We analyze the consequences of different hydro- 
dynamic resolutions on the nucleosynthesis, we compare their 
different distribution of burned and unburned material, we dis- 
cuss the trend of the Y^, in the markers as a consequence of our 
nucleosynthesis calculations. Finally, we discuss the velocity 
distribution of different nuclear species. 
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4.1. Nucleosynthesis and yields: comparison between 
2D and 3D 

For the nucleosynthesis calculations the peak temperatures 
combined with the density distribution achieved during the 
propagation of the front, are the most important quantities. As 
shown in Fig. 5, for the b303d_768 model the maximum of 
the T in the markers covers a large range (1.5x10^ K< T < 
8.4x 10^ K). Note that in this Figure each dotted line represents 
a tracer particle trend (we selected randomly some markers for 
this plot), and the thick line is the upper envelope of the T and 
p distribution in the total tracer particle sample. At high tem- 
peratures all strong interactions and photodisintegrations are 
so fast that a chemical equilibrium (nuclear statistical equilib- 
rium, NSE) is immediately achieved (in our calculations we as- 
sume NSE condition for T > 6.0xlO^K). The resulting chem- 
ical composition is therefore only dependent on p, T, and the 
neutron eccess (that is determined by the total amount of elec- 
tron captures taking place on longer timescales). An example 
of the behaviour of the chemical abundances in one marker as 
a function of time is shown in Fig. 6. For this case the marker is 
originally located at a radius of ~ 150 km, i.e. in the innermost 
dense zone; the initial density at the position of the marker is ~ 
2.5x10^ gr/cm-' (the central density for this model is 2.9x10"^ 
gr/cm^) and it is heated by the flame front almost immedi- 
ately (at ~0.1 s). The NSE conditions are achieved very fast 
and the temperature reaches a quite high peak T ~8.5xlO'^K 
(Fig. 6, upper left panel). The resulting electron fraction Y^, 
drops rapidly (in ~0.3 s) from the initial value of 0.4989 to 
0.468 (Fig. 6, upper right panel). The main O-burning prod- 



ucts. Si and S are shown in Fig. 6 (left middle panel), to- 
gether with the abundance of ^"^Mg resulting from C-burning. 
Furthermore, the rise of the temperature in excess of 6x10^ 
K leads to a complete NSE and ^^Si exhaustion within 0.1 s. 
Due to the high density (typical of explosive Si-burning) nor- 
mal freeze-out occurs. One of the main products is ^*Fe (Fig. 6, 
right middle panel). Since T and p conditions are very high, 
also neutron-rich nuclei are built up due to electron captures, 
and ^^Fe is partly replaced by ^"^Fe (Fig. 6, right middle panel) 
and ^^Ni (Fig. 6, left lower panel). In such details we can fol- 
low the nucleosynthesis changes along the time in all the tracer 
particles. 

In Fig. 7 we show the distribution of Yj, vs. time obtained 
as a result of the nucleosynthesis calculation in the model 
b30^dJ768. The two thick lines stand for the upper and lower 
values of the Y^, in the markers. The dashed lines represent the 
Y(, time evolution in some markers randomly taken as exam- 
ples. The range covers values from ~0.5 (that represent the typ- 
ical initial composition as well as the composition of markers 
with rather low T), down to ~0.462, reached by the markers 
with the highest T (see e.g. the example in Fig. 6). 

Table 1 lists the synthesized masses for all the stable 
isotopes up to ^^Zn for the models c3J2d^l2, c3^dJ256, 
b5_3dJ256, and b30_3dJ768. For comparison we also include in 
column 2 the calculations for the W7 model (from Thielemann 
et al. 2003). For the calculations presented in column 5,6,7 
we include the nucleosynthesis results starting only when the 
temperature has reached 90% of the peak. In fact, due to fi- 




0.6 1 

tinie(s) 



Fig. 6. Example of the nucleosynthesis calculation in one tracer 
particle. T and p are plotted in the upper left panel; with T > 
6x10*^ K reaches NSE conditions. In the upper right panel the 
resulting Y^ is shown. The other panels give the mass fraction 
vs. time for ""O, ^^Si, ^^S, ^'^■^^Fe, 56-57,58j^j^ 63^^^^ ^j^^j 642ji. 
For '^O, ^'^Si, ^^S the time is plotted only up to 0.2 s since their 
abundances are zero at later times. 



nite numerical resolution in the hydrodynamic simulation, the 
rise of the temperature is not as steep as it would be in re- 
ality. Consequently, when markers reach NSE conditions and 
weak-interactions start to play the most important role, the nu- 
cleosynthesis timescales are fast, and even 0.1 s (i.e. typical 
timescale we have for the rise of T) are crucial for some re- 
actions to give an important contribution. Nevertheless, as one 
can see from the Table comparing column 4 (the nucleosyn- 
thesis for the c3^dJ256 model has been calculated considering 
the all rise of the T curve) and column 5 (the nucleosynthesis 
for the same c3^dJ256 model has been calculated only when 
T has reached 90% of the peak), the differences in the total 
yields are very small. This is due to the fact that the amount 
of markers aff'ected by this inaccuracy is a small fraction of 
the total. In the Table we do not include isotopes heavier than 
^^Zn, even if the network we used was extended up to ^^Mo, 
since their resulting mass fraction are smaller than ~1O"'^M0. 
An important thing to notice is the difference in the amount 
of unburned material (as defined at the beginning of this pa- 
per, we consider unburned the material with T < 1.5x10'' K) 
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Fig. 7. The same as Fig. 5 for the history of the Y^. in the mark- 
ers. 

in the four models. For the 2D model c3_2d-512 ~60% of 
the total material remains unburned, instead of 40% we ob- 
tain for the 3D model c3_3d-256. Both of these models are 
centrally ignited, therefore the difference in the amount of un- 
processed material is mainly a consequence of the difference 
in the total energy distribution due to multi-dimension effects 
(see Fig. 1). In Fig. 8 we show the distribution of burned and 
unburned material for the model b30^dJ768 (upper panel), 
c3^d-256 (middle panel), and b5^d-256 (lower panel), in 
terms of number of tracer particles. The distribution plotted 
in Fig. 8 has been taken at 1.1 sec for all the three models, 
that is the final time we reached in the b30^dJ768. For all 
these three simulations the burned component dominates at a 
radius of ~5.0xl0^cm, instead the unburned material is more 
or less uniformly distributed, with a dominant component in 
the outermost zones and a tail in the central regions. As also 
discussed by Reinecke et al. (2002b) the distribution of the 
unprocessed material depends on the initial conditions for the 
burning. When the model is centrally ignited, as c3^d-256, 
the dominant component of unburned material is in the outer- 
most regions. Instead with a floating-bubble model with com- 
parable resolution, like b5Jd^56, unprocessed material can 
also be concentrated in the center Nevertheless the amount of 
unprocessed material in the center for a floating-bubble model 
depends on the amount of ignition-spots together with the reso- 
lution used. In fact (see Fig. 8) for the b30^dJ768 case most of 
the unburned component in the innermost regions disappeared. 
We also tested consequences of burning for longer times. 
As mentioned in Section 3.2, at ~1.2 sec the upper limit of 
the temperature is ~4xl0^ K (see also Fig. 5), and explosive 
C- and a bit of Ne-burning can still occur. Therefore, to calcu- 
late the nucleosynthesis we extrapolated temperature and den- 
sity for ~1 sec. further (when the upper values for T are not 
higher that 1.5x10^ K). If we compare the nucleosynthetic 
yields calculated at the end of the hydro simulation with the 
nucleosynthetic yields calculated using the extrapolated val- 
ues of T and p, we obtain that only few % of '^C burns. The 
consequences can be relevant only for few isotopes, the main 
products of explosive C-burning, i.e. ^"'^'Ne, ^^Na, ^^•^''Mg, 



b30_3d_768 - 




2.0 108 4.0 108 6.0 lO^ 

R [cm] 

Fig. 8. Distribution of the tracer particles vs. radius at ~ 1 .2 sec. 
for b30^dJ768 (upper panel), c3^d-256 (middle panel), and 
b5^dJ256 (lower panel) models. The unburned particles (T < 
1.5x10^ K) are plotted with a solid line, and the processed par- 
ticles with a dashed line. 



^^Al. Therefore only for those isotopes we expect important 
changes if we could follow the hydrodynamic simulation for 
longer timescales. 

It is interesting to notice the difference in the '^C and '^O 
abundance of the W7 model and our multi-dimensional mod- 
els. The '^C of the W7 is about a factor of ~10 lower than 
in our cases, instead "'O of W7 is only a factor of ~3 lower. 
While "'O and '^C we obtain are built up by unprocessed par- 
ticles, in the W7 C-burning is more efficient and burns a signif- 
icant amount of '^C at low T and p, with a resulting different 
C/O ratio. Also, the ^*'Fe mass (mainly deriving from the decay 
of the long-Uved ^''Ni) obtained in the W7 model (0.696 Mq, 
Thielemann et al. 2003) is by far higher than the ^^Fe mass re- 
sulting from the multi-dimensional SNIa models. The highest 
value we can reach in our models is obtained with the highest 
resolution floating-bubble model b30JdJ68 (0.44 Mq). We 
note that the initial conditions (in this case 30 ignition spots) 
are crucial for a more precise study of the nucleosynthesis, in 
particular of the innermost regions. Possibly, a model with even 
higher number of ignition bubbles would give us a still higher 
^''Fe mass. 

We note that neutron-rich isotopes like "^^Ca, ^"Ti, ^"^Cr, 
54pg^ 58|^j ^g strongly underproduced with respect to the W7 
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Fig. 9. Nucleosynthetic yields (in mass fraction normalized to the solar value and to the corresponding solar ratio) obtained using 
19683 tracer particles in the 3D model b30^dJ768 compared to the W7 yields given by Thielemann et al. (2003). 



model presented by Iwamoto et al. (1999), instead are in good 
agreement with the Branchwitz et al. (2000) and Thielemann 
et al. (2003) predictions. This is due to the differences in the 
electron capture rates adopted (as just discussed above). As a 
consequence, electron capture rates of nuclei affects directly 
the electron fraction Y,, (Thielemann, Nomoto, & Yokoi 1986). 
As shown in Fig. 7, the lowest Yg reached in the markers is 
0.462, instead the lowest value in the W7 model by Iwamoto et 
al. (1999) is ~ 0.446. For the model by Iwamoto et al. (1999) 
low Ye values (< 0.46) are reached in the innermost zones, 
i.e. M < 0.03Mq, where the highest temperatures are reached 
(~9.0xl0^ K). The highest T reached in our models, as shown 
in Fig. 5, are ~8.4xl0^ K with consequently higher Yg. 

Finally, in Table 2 we report the synthesized masses of the 
main radioactive species from ^^Na up to ''•'Ni. Bigger differ- 
ences between our models and W7 are the abundances of "^^Ca 
and ^*'Zn (in this case we use for comparison the W7 data from 
iwamoto et al. 1999, since the Thielemann et al. (2003) data 
are not yet available for the long-lived isotopes). As discussed 
in Thielemann et al. (2003), these isotopes are very sensitive to 
small variations in the central density of the model. We used 
2.9x10^ gr/cm"* instead of the 2.0x10^ gr/cm^ of the model 
B2C20 presented by Thielemann et al. (2003) and used in the 
current paper as a comparison to our model. 

In Fig. 9 we show the yields obtained for our 'standard' 
model b30JdJ768 compared with the W7 yields (Brachwitz et 
al. 2000, Thielemann et al. 2003), scaled to their relative solar 
abundances and to the ^''Fe abundance. As one can see from 
the figure, a part from the difference in the relative ^*'Fe mass 
(0.44 Me for the b30JdJ68 and 0.669 Mq for W7), as well as 



in the unburned material (i.e. '^C, '*0 and ^^Ne), the trend for 
the production of different isotopes is quite similar. 

In Fig. 10 we plotted the yields of the c3J2d-512 and 
c33dJ256 {upper panel), and b5^d-256 (lower panel), nor- 
malized to the 'standard' model b30^dJ768. From Fig. 1 1 one 
can notice that with a similar amount of burned material (the 
difference between the two models in the total burned mate- 
rial it is not more than 10%), the b5^d-256 produces much 
more a-elements. Infact the higher efficiency of '^C burning is 
clearly shown by a higher production of ^^Na and "^''Ca. Also 
the b5^d-256 model has a lower ^^Fe production (0.34 Mq, in- 
stead of 0.44 Mq of the b30^dJ768), and Fe-group elements. 

4.2. Radial velocity distribution 

In Fig. 1 1 we show the mass fractions of selected isotopes as a 
function of the radial velocity, taken at the end of our simula- 
tion (i.e. ~1.2 sec.) for the model b30^dJ768 and c3_3d-256. 
With thick lines we plot the unprocessed material in the form 
of '^C, '^O, and ^^Ne. We note that at the lowest velocities 
(~1000 km/s) the dominant component is represented by the 
unburned material (i.e. '^C and '^O) for the floating-bubble 
high-resolution model, and by ^^Fe for the centrally ignited 
model. On the opposite site, at the highest velocities (> 10000 
km/s) the unburned material dominates in the centrally ig- 
nited model, instead is comparable to the ^^Fe component in 
the floating-bubble high resolution model. Maximum veloci- 
ties reached are about 12000 km/s in both cases. As mentioned 
in the previous Section, our models at 1.2 sec are not yet in ho- 
mologous expansion, i.e. pressure and gravity still play a role. 
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Fig. 10. Nucleosynthetic yield ratio, comparing the model 
c3J2.d-512 (dotted line), and c3_3d-256 (solid line), with our 
standard model b30_3dJ768 (upper panel). Nucleosynthetic 
yield ratio, comparing the model b5^dJ256 with our standard 
model b30_3dJ768 (lower panel). 

changing the velocity distribution and possibly also the den- 
sity, therefore the velocity distribution of the elements has to 
be taken with care. However the distribution in velocitites 
might leave observable features in the spectra that could be 
used for diagnostic purposes. Finally, by projecting our 3D 
results on previously obtained ID ones, the observed evolu- 
tion of the Si, S, and Ca Unes might give us the possibiUty 
to understand whether normal SNIa are well mixed defla- 
grations. 

5. Summary and conclusions 

In this paper we presented the results of nucleosynthesis cal- 
culations obtained coupling a tracer particle method to two- 
dimensional and three-dimensional Eulerian hydrodynamic 
calculations of SNIa explosion. The multidimensional SN la 
simulations described in this work employed a minimal nuclear 
reaction network, sufficient for a good approximation of the 
thermonuclear energy release. Although the predicted chemical 
composition agrees well with the expectations, we presented 



here the results of a very detailed study of the nuclear abun- 
dances in the ejecta obtained by post-processing the output 
of the hydrodynamic models. Since the multidimensional hy- 
drodynamical scheme applied is of Eulerian type (i.e. the grid 
does not move with the fluid), we added a lagrangian compo- 
nent to the calculations in the form of tracer particles. In order 
to record temperature and density evolution as a function of 
time (necessary input for the nucleosynthesis calculations) we 
homogeneously distributed ~20000 marker particles and fol- 
lowed their T and p evolution. We then calculated the nucle- 
osynthesis experienced by each marker and computed the total 
yield as a sum over all the markers including the decays of un- 
stable isotopes. 

The nuclear reaction network employed in computing the 
explosive nucleosynthesis contains 383 nuclear species rang- 
ing from neutrons, protons, and a-particles to molibdenum. 
For this work, the initial mixture we used consists of (mass 
fraction) 0.475 Mq of '^c, 0.5 Mq of '''O, and 0.025 Mq of 
^^Ne. When the flame passes through the fuel, C, O and Ne 
are converted to heavier elements, with different compositions 
depending on the T and p history. Nuclear statistical equilib- 
rium conditions are assumed in the marker particles with T > 
6 10^ K. At such temperatures (r > 6 10^ K) a mixture of ^''Ni 
and a-particles in NSE is synthesized. Below that temperature 
burning only produces intermediate mass elements. Once the 
temperature drops T < 1.5 10^ K, no burning takes place dur- 
ing the short timescale (^ 1.5 s) of the explosion ("unburned" 
material). 

The current research focused on the sensitivity of the ex- 
plosion on the ignition conditions and on the detailed nucle- 
osynthetic yields that they predict. We could demonstrate that 
multi-dimensional explosion models allow us to predict their 
nucleosynthesis yields with some confidence. It was shown that 
only 3D models are potentially able to produce enough ^''Ni to 
explain the light curves of "normal" type la supernovae, and 
that also the ignition conditions (central ignition vs. several off- 
center ignition spot) affect the nucleosynthesis yields. Since the 
number of ignition spots we can put into the numerical models 
depend on the spatial resolution and since the explosion en- 
ergy as well as the Ni-mass increase with increasing number 
of spots, we expect that our best resolved b30^dJ768 model is 
closest to what a "typical" pure-deflagration supernova might 
eject. The general nucleosynthesis outcome of SNela is dom- 
inated by Fe-group elements, involving also sizable fractions 
of Si-Ca and minor amounts of unburned (C and O) or pure 
C-burning products (e.g. Na, Ne, Mg). Despite of the fact that 
differences with respect to the standard W7 nucleosynthesis 
(Iwamoto et al. 1999, Brachwitz et al. 2000, Thielemann et 
al. 2003) are found, in particular in the ^^Ni mass produced, as 
well as in the final amount of unburned material, in general the 
nuclear yields are consistent with expectations. We can there- 
fore say with some confidence and without parameterization, 
that the Chandrasekhar mass scenario with a pure turbulent de- 
flagration is a viable candidate for SN la explosions. We also 
note that the significant amount of unburned material ejected 
by our SNIa models may have an interesting impact on the role 
of SNIa in the context of Galactic chemical evolution of C (giv- 
ing a contribution of the order of ~20% to the total C at the 
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solar composition). In the case of O still the main sources are 
massive stars. 

Comparing the nucleosynthesis presented in this paper to 
observed SNIa spectra, the reader should keep in mind that our 
models do not reach the homologous expansion phase. We are 
currently working to modify the combustion hydrocode, us- 
ing a moving grid that will allow us to follow the evolution 
much longer. The results will be published elsewhere. We are 
also performing a detailed parameter study of the variation of 
the central density and of the initial carbon/oxygen ratio of the 
SNIa progenitor (Ropke et al., in preparation). Finally, recent 
calculations by Timmes et al. (2003) indicate large variations 
of the ^^Ni mass as a function of metrallicity (measured by the 
original ^^Ne content). An investigation of the metallicity effect 
on the nucleosynthesis and yields is also in progress (Travaglio 
et al., in preparation). 
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Fig. 11. Mass fractions of selected isotopes as a function of the radial velocity of the markers (taken at ~1.2 sec) for the 
b303d-768 model (upper panel) and for the c3_3d-256 (lower panel). The width of each velocity bin is 300 km/s. For each 
isotope we sum its abundance over all markers in a certain velocity bin. 
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TABLE I 

Synthesized Mass (M©) in SNIa models 



Species 


W7<-->' 


cSJdjW-"' 


c3JdJ5ff''> 


c3JdJ56^'> 


b5JdJ56'-'> 


b30JdJ6»'> 


l^C 


5.04E-02 


4.09E-01 


3.37E-01 


3.37E-01 


2.49E-01 


2.78E-01 


i^C 


1.07E-06 


1.07E-10 


9.71E-11 


9.68E-11 


8.21E-06 


3.98E-06 


"N 


4.94E-07 


2.71E-09 


3.96E-09 


3.48E-09 


1.04E-03 


2.76E-04 


15n 


1.25E-09 


4.40E-11 


7.18E-11 


6.99E-11 


2.48E-05 


1.23E-06 


160 


1.40E-01 


4.74E-01 


4.16E-01 


4.17E-01 


3.90E-01 


3.39E-01 


"0 


3.05E-08 


1.16E-09 


1.29E-09 


1.13E-09 


7.81E-06 


1.31E-06 


180 


7.25E-10 


9.49E-11 


1.62E-10 


1.52E-10 


1.15E-04 


l.OlE-05 


Wp 


5.72E-10 


2.64E-11 


3.72E-11 


3.34E-11 


1.08E-06 


2.84E-08 


2»Ne 


1.97E-03 


4.70E-03 


7.39E-03 


7.10E-03 


3.18E-02 


6.28E-03 


2'Ne 


8.51E-06 


7.11E-07 


1.14E-06 


1.03E-06 


5.96E-05 


2.16E-05 


22Ne 


2.27E-03 


2.15E-02 


1.77E-02 


1.77E-02 


1.14E-02 


1.42E-02 


2^Na 


6.20E-05 


2.99E-05 


5.09E-05 


5.10E-05 


3.49E-03 


8.65E-04 


^"Mg 


1.31E-02 


1.04E-02 


1.48E-02 


1.26E-02 


2.35E-02 


7.53E-03 


^'Mg 


4.71E-05 


5.49E-05 


8.57E-05 


7.64E-05 


2.41E-03 


5.13E-04 


26Mg 


3.31E-05 


6.60E-05 


1.06E-04 


l.OlE-04 


8.56E-04 


1.81E-04 


"Al 


8.17E-04 


7.39E-04 


1.08E-03 


9.73E-04 


2.11E-03 


5.85E-04 


2*Si 


1.52E-01 


4.42E-02 


5.89E-02 


5.39E-02 


1.19E-01 


5.39E-02 


^'Si 


7.97E-04 


6.47E-04 


9.49E-04 


9 27F-04 


1.81E-03 


5.61E-04 


3»Si 


1.43E-03 


1.06E-03 


1.48E-03 


1.31E-03 


2.20E-03 


8.03E-04 


31p 


3.15E-04 


2.02E-04 


2.85E-04 


2.69E-04 


5.24E-04 


1.72E-04 


32 s 


8.45E-02 


1.60E-02 


2.22E-02 


2.57E-02 


5.70E-02 


2.62E-02 


33s 


4.11E-04 


1.05E-04 


1.42E-04 


1.58E-04 


3.21E-04 


1.21E-04 


34s 


1.72E-03 


8.68E-04 


1.15E-03 


1.15E-03 


2.30E-03 


1.04E-03 


36 s 


2.86E-07 


1.64E-07 


2.24E-07 


2.47E-07 


3.95E-07 


1.53E-07 


35C1 


1.26E-04 


3.60E-05 


4.88E-05 


5.90E-05 


1.31E-04 


4.58E-05 


3'Cl 


3.61E-05 


6.89E-06 


8.97E-06 


1.27E-05 


3.21E-05 


1.21E-05 


^OAr 


1.49E-02 


2.12E-03 


3.14E-03 


4.09E-03 


9.04E-03 


4.24E-03 


3«Ar 


8.37E-04 


3.30E-04 


4.13E-04 


5.12E-04 


1.20E-03 


5.59E-04 


""Ar 


1.38E-08 


1.49E-09 


2.06E-09 


3.04E-09 


4.92E-09 


1.91E-09 


39k 


6.81E-05 


1.51E-05 


1.84E-05 


2.95E-05 


7.69E-05 


3.24E-05 


«K 


6.03E-06 


9.03E-07 


1.17E-06 


2.20E-06 


6.03E-06 


2.41E-06 


""Ca 


1.21E-02 


1.68E-03 


2.66E-03 


3.40E-03 


7.08E-03 


3.59E-03 


«Ca 


2.48E-05 


6.66E-06 


8.43E-06 


1.41E-05 


3.61E-05 


1.58E-05 


«Ca 


1.07E-07 


2.26E-08 


3.06E-08 


3.96E-08 


6.37E-08 


5.10E-08 


«Ca 


9.62E-06 


1.80E-06 


2.81E-06 


3.10E-06 


4.52E-06 


3.61E-06 


""Ca 


2.44E-09 


2.58E-12 


3.46E-12 


1.14E-11 


1.91E-11 


8.53E-12 


«Ca 


1.21E-12 


1.99E-17 


3.20E-I7 


1.05E-16 


1.54E-16 


4.01E-15 


«Sc 


2.17E-07 


2.16E-08 


3.06E-08 


6.08E-08 


1.65E-07 


6.47E-08 


«Ti 


1.16E-05 


2.80E-06 


3.53E-06 


5.62E-06 


1.47E-05 


6.62E-06 


47.J.J 


5.45E-07 


1.38E-07 


1.88E-07 


2.20E-07 


3.61E-07 


2.64E-07 


48Ti 


2.07E-04 


4.11E-05 


6.96E-05 


7.28E-05 


1.32E-04 


7.69E-05 


«Ti 


1.59E-05 


3.28E-06 


5.22E-06 


5.62E-06 


l.lOE-05 


5.78E-06 


'"Ti 


1.62E-06 


8.22E-10 


2.08E-08 


2.08E-08 


8.67E-10 


2.67E-07 


50 Y 


4.58E-09 


2.04E-09 


3.66E-09 


3.50E-09 


2.69E-08 


2.66E^09 


51v 


3.95E-05 


1.70E-05 


1.90E-05 


1.90E-05 


2.89E-05 


1.95E-05 


'"Cr 


2.23E-04 


1.20E-04 


1.20E-04 


l.lOE-04 


1.67E-04 


1.19E-04 


52Cr 


4.52E-03 


1.91E-03 


2.80E-03 


2.76E-03 


3.48E-03 


2.58E-03 


53 Cr 


6.49E-04 


4.78E-04 


5.18E-04 


4.81E-04 


5.09E-04 


4.83E-04 


'"Cr 


3.04E-05 


3.42E-06 


6.38E-06 


5.92E-06 


4.11E-06 


1.22E-05 


^^Mn 


6.54E-03 


5.63E-03 


5.93E-03 


5.53E-03 


4.53E-03 


6.38E-03 


^■•Fe 


7.49E-02 


6.79E-02 


6.61E-02 


6.21E-02 


4.48E-02 


7.33E-02 


^-ipe 


6.69E-01 


2.44E-01 


3.28E-01 


3.36E-01 


3.40E-01 


4.39E-01 


"Fe 


2.52E-02 


1.05E-02 


1.35E-02 


1.36E-02 


1.28E-02 


1.86E-02 


5«Fe 


1.74E-04 


8.25E-06 


3.16E-05 


3.02E-05 


8.58E-06 


1.05E-04 


^'Co 


7.66E-04 


6.70E-04 


7.62E-04 


6.81E-04 


4.53E-04 


7.33E-04 


5«Ni 


1.02E-01 


6.13E-02 


7.52E-02 


7.31E-02 


5.56E-02 


9.66E-02 


""Ni 


9.22E-03 


7.23E-03 


9.24E-03 


8.16E-03 


5.39E-03 


7.73E-03 


«Ni 


2.69E-04 


6.11E-05 


8.86E-05 


9.26E-05 


9.99E-05 


1.13E-04 


62Ni 


2.31E-03 


5.71E-04 


7.78E-04 


8.16E-04 


9.21E-04 


1.12E-03 


"Ni 


1.84E-07 


2.73E-11 


1.61E-09 


1.61E-09 


1.93E-10 


5.29E-08 


"Cu 


I.59E-06 


9.24E-07 


9.26E-07 


9.27E-07 


8.20E-07 


9.56E-07 


O'Cu 


7.72E-07 


1.88E-07 


2.51E-07 


2.61E-07 


2.81E-07 


3.77E^07 


"Zn 


1.50E-05 


3.72E-06 


4.47E-06 


4.65E-06 


4.83E-06 


6.78E-06 


«'Zn 


1.31E-08 


6.11E-06 


7.55E-06 


7.86E-06 


8.90E-06 


1.16E-05 


"Zn 


1.18E-11 


4.15E-09 


5.49E-09 


5.68E-09 


6.53E-09 


7.96E09 


""Zn 


2.66E-10 


2.85E-09 


3.68E-09 


3.86E-09 


4.86E-09 


5.26E-09 



Thielemann et al. (2003) 

This work 

This work. In this run we allow the nucleosynthesis calculations for those tracer particles that reach NSE conditions only starting at 90% of the temperature peak. 
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TABLE 2 
Synthesized Mass (Mq) for radioactive species in SNIa models 



Species 


W7<^' 


cs.icun'-'''' 


c3JdJ.5&^^ 


c5.id^J6<'-) 


b5JdJ.56^''^ 


biojdje^''^ 


^^Na 


1.73E-08 


5.76E-08 


9.42E-08 


8.94E-08 


8.00E-07 


l.OOE-07 


2«A1 


4.93E-07 


5.98E-07 


9.82E-07 


9.11E-07 


4.53E-05 


4.47E-06 


3«C1 


2.58E-06 


5.56E-07 


7.62E-07 


9.73E-07 


1.99E-06 


6.32E-07 


s'Ar 


1.20E-08 


1.68E-09 


2.34E-09 


3.24E-09 


9.41E-09 


2.24E-09 


40k 


8.44E-08 


7.77E-09 


l.lOE-08 


1.73E-08 


4.10E-08 


1.23E-08 


^'Ca 


6.09E-06 


9.01E-07 


1.17E-06 


2.20E-06 


6.02E-06 


2.40E-06 


44-pi 


7.94E-06 


1.80E-06 


2.80E-06 


3.09E-06 


4.50E-06 


3.61E-06 


48 Y 


4.95E-08 


4.10E-05 


6.95E-05 


7.27E-05 


1.32E-04 


7.68E-05 


49y 


1.52E-07 


3.28E-06 


5.22E-06 


5.62E-06 


l.lOE-05 


5.78E-06 


"Mn 


2.77E-04 


4.77E-04 


5.16E-04 


4.78E-04 


5.09E-04 


4.79E-04 


60pg 


7.52E-07 


1.44E-14 


6.57E-12 


6.57E-12 


8.54E-13 


4.52E-10 


«Co 


1.44E-04 


9.50E-05 


1.04E-04 


9.81E-05 


1.69E-04 


1.32E-04 


"Co 


1.48E-03 


1.33E-03 


1.38E-03 


1.16E-03 


6.35E-04 


1.15E-03 


«»Co 


4.22E-07 


4.00E-10 


6.82E-09 


6.77E-09 


1.12E-09 


2.66E-08 


«Ni 


5.86E-01 


2.16E-01 


2.95E-01 


3.08E-01 


3.27E-01 


4.18E-01 


"Ni 


2.27E-02 


9.17E-03 


1.21E-02 


1.25E-02 


1.22E-02 


1.74E-02 


'"Ni 


6.71E-04 


6.62E-04 


7.40E-04 


6.60E-04 


4.47E-04 


7.11E-04 


63Ni 


8.00E-07 


4.15E-11 


2.10E-09 


2.10E-09 


1.98E-10 


2.22E-08 



" Iwamotoetal. (1999) 
'' This work 

'^ This work. In this run we allow the nucleosynthesis calculations for those tracer particles that reach NSE conditions only starting at 90% 
of the temperature peak. 



